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The magnetic anisotropy energy of Fe and Co monatomic nanowires both free-standing and at 
the step edge of a Pt surface is investigated within the framework of the density-functional theory 
and local-spin-density (LSDA) approximation. Various types of the analysis of the calculated data 
reveal that the spin-orbit interaction of the Pt atoms and the hybridization between the electronic 
states have an important impact on the direction of the easy axis and on the magnitude of the 
magnetic anisotropy, both by a direct contribution localized at the Pt atoms and by an indirect 
contribution due to the modification of the Co-localized part via hybridization effects. 
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I. INTRODUCTION 



The experimentally observed existence of ferromag- 
netism in linear chains of Co atoms, grown at a step 
edge of a Pt(997) surface terracei, motivated several 
theoretical investigations which were aimed to interpret 
the measured large orbital contribution to the magnetic 
moments 2 ^, and the magnetic anisotropy 4 ' 5 - 6 7 - 8 - 9 . The 
latter is distinguished by the experimentally found easy 
axis along the direction perpendicular to the Co wire, 
shifted by 43° from the surface normal towards the Pt 
step edge, and by the magnetic-anisotropy energy (MAE) 
of the order of 2 meV/Co atom. This value is a fac- 
tor of about 50 larger than the one of hep Co (which 
is already large for a transitional metal). Various theo- 
retical attempts to reproduce the experimental findings 
yielded quantitatively different results. Hong and Wu^ 
performed the full-potential linearized-augmented-plane- 
waves (FLAPW) calculation for infinite Co wires on the 
flat Pt(001) and Cu(001) substrates. They found that 
the MAE magnitude was enhanced when they replaced 
the Cu substrate by platinum due to the strong spin- 
orbit coupling (SOC) of Pt atoms. They obtained the 
easy axis along the wire in the case of the Cu(001) sub- 
strate, and perpendicular to the wire in the plane of the 
Pt(001) substrate. They also noticed a certain sensitiv- 
ity of the results on the size of the supercell. The results 
of the Korringa-Kohn-Rostoker (KKR) calculations 5 - 6 on 
finite (up to 10 atoms) Co wires along a step edge of a 
Pt(lll) surface were in remarkable agreement with the 
experiment. However, to some extent, this agreement 
might be coincidental since the atomic-sphere approxi- 
mation (ASA) was applied for the potential. Again, the 
authors ascribed the peculiar magnetic anisotropy to the 
influence of the Pt substrate. The same conclusion was 
drawn also on the basis of the FLAPW calculations for 
the infinite Co 7 * 8 * 9 - and Fem wires at the Pt(lll) step edge 
by analyzing the data in the frame of a symmetry-based 
"atomic-pair" model. Although shape approximations 
were not applied for the potential, the experimental easy 



axis was not reproduced as good as in the case of the 
KKR ASA calculation 5 * 6 -, even if the structural relax- 
ations were taken into account^. Finally, an FLAPW 
calculation^ for a considerably larger supercell as used 
m Ref. H including structural relaxations, found an easy 
axis which is in a plane parallel to the substrate sur- 
face and almost perpendicular to the wire, in contrast to 
the result of Ref. and to the experimental result. All 
this indicates that the orientation of the calculated easy 
axis depends extremely sensitively on the calculational 
details. We therefore do not take the enormous effort 
to converge the results with respect to all parameters of 
the calculation (e.g., geometry and size of the supercell), 
but we concentrate on the physical interpretation of the 
results concerning the influence of the Pt substrate by 
means of two different types of the analysis. 



II. CALCULATIONAL METHOD 

We performed the calculations with the FLAPW 11 
WienQTi 2 - code by using the local-spin-density approxi- 
mation (LSDA) 1 ^ for the exchange-correlation potential. 
The SOC contribution was added to the Hamiltonian 
in terms of the second-variational method 1 ^ 5 -. In the 
implementation of Ref. 12 SOC is taken into account 
just within muffin-tin spheres but not in the interstitial 
region. The energy difference in the total energy due 
to different magnetization directions was determined ac- 
cording to the magnetic force theoremi&ii. In this ap- 
proximation, first a self-consistent effective potential 
is determined for an initial magnetic configuration (0) of 
the system. In a subsequent step the eigenvalues are de- 
termined for respectively two different orientations (1) 
and (2) of the magnetization by respectively one single 
diagonalization of the Hamiltonian for the same fixed po- 
tential V^g' but with the SOC term corresponding to 
the two orientations of the magnetization. The magnetic 
anisotropy energy eMCA which is defined as the difference 
in the total energy for the two magnetization directions 
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is then approximated by the difference in the respective 
sums of eigenvalues: 



e MCA 



eM (2) (e) de - 



eM {1) {e)de, (1) 



In eq. (1) the functions Af^ (e) and (e) are the 
electronic density of states for the two orientations, and 



and E%> are the respective Fermi levels, which are 
calculated from the number Z of the valence electrons 
according to: 



p(2) 



z = 



(e) de. 



(2) 



The convergence of the magnetic anisotropy energy with 
respect to the number of k points used for the sampling of 
the Brillouin zone is notoriously bad. However, it turns 
out that less k points are needed for the calculation of 
v}2\ Hence, the advantage of the magnetic force the- 
orem is that the calculation of e^ CA from (1) requires 
considerably less computational power than the deter- 
mination of the respective difference in the total energy, 
obtained from two self-consistent calculations. Because 
the FLAPW method is very time consuming, the applica- 
tion of the force theorem is desirable. An open question 
thereby is how to select the initial configuration (0) for 
the determination of V^g in order to minimize the dif- 
ference between cmca and e^ CA . To figure this out we 
have calculated by means of the linear-muffin-tin-orbital 
(LMTO) method in atomic-sphere approximation 18 both 
eMCA (which is possible because LMTO-ASA is much 
faster, albeit less accurate, than FLAPW) and e^ CA for 
various configurations (0). It turns out that it is abso- 
lutely indispensable to include SOC for the initial con- 
figuration. The initial orientation of the magnetization 
has only a subtle effect on e^ CA since for each orienta- 
tion the agreement between eMCA and e^ CA is nearly 
perfect as presented in Fig. 1. For this test, the infinite 
Co nanowire at the step edge of the Pt(lll) surface was 
modelled by a supercell with the geometry from Ref. 0, 
which contains 16 atoms (15 Pt + one Co atom) and 8 
empty atomic spheres that are required due to the ASA 
approximation. 
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FIG. 2: The calculated energy difference e~MCA = e ( 
0°,(p) - e(8 = 0°,(p= 0°) for the Co (top) and the Fe (bot- 
tom) Pt-supported wire with SOC at all sites + or just at the 
Co(Fe) x or Pt □ atoms. 



The rest of the calculations were performed for a 
monatomic wire at the step edge of a Pt surface described 
by a supercell of 13 atoms (12 Pt + one Co or Fe atom) 
by adopting the experimental lattice parameters of fee 
Pt (for details, see Ref. 0). In addition, freestanding 
wires were considered by removing all Pt atoms from the 
supercell. The orientation of the coordinate system was 
adopted from Refs 7 8 9 (see the inset in Fig. 2). The 
cut-off parameter for the plane- wave expansion was set to 
7.3 Ry. The Brillouin-zone (BZ) integration was carried 
out by both the modified tetrahedron— and the Gaus- 
sian smearing 20 methods with 45 (40 in the case of the 
freestanding wires) k-points in the full BZ for the self- 
consistent part and 192 (184) for the force-theorem part 
of the calculation as determined on the basis of conver- 
gence tests. 

We calculated e^cA f° r the magnetization confined ei- 
ther to the yz (the angle 0, Fig. 2) or the xz plane (the 
angle 8, Fig. 3) with respect to the energy for the case 
of the magnetization along the z axis. The calculated 
data were fitted by the functions a — b cos 2 (<fi — c) and 
a — fecos 2 (8) for the respective planes. Since there is a 
quantitative deviation from the results obtained by Shick 
et ali& for a slightly different supercell, we performed a 
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FIG. 3: The calculated energy difference Cmca = e (#i 4 1 — 
0°) - e(6 = O°,0 = 0°) for the Co (top) and the Fe (bot- 
tom) Pt-supported wire with SOC at all sites + or just at 
the Co(Fe) x or Pt □ atoms. The symbols ■ present the 
respective energy difference for the freestanding wires. 
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FIG. 4: The calculated magnetic anisotropy energy e£ AE = 
e(6 = 0°,(p = 22.85°) - e{6 = 0°,(j> = -67.15°) for the Co 
wire (top) and e M AE = e(9 = 0°, <j> = 80.40°) - e(0 = 0°, <j> = 
—9.60°) for the Fe wire (bottom); total contribution (+), Co 
(Fe) contribution (x), Pt contribution (□). 



test calculation with their type of the supercell (16 atoms 
but considerably less vacuum at the top of the (111) sur- 
face) , and we could satisfactorily reproduce the published 
results from Refs. L^la- This demonstrates that different 
computer codes give the same results, and it confirms the 
findings from previous investigations that the magnetic 
anisotropy depends extremely sensitively on the details 
of the structural model. 

To explore the influence of the Pt substrate we per- 
formed two types of the analysis: 

(a) The densities of states AA 1 ' 2 -* (e) in eq. (1) can be 
subdivided into the contributions A/1 1 ' 25 (e) which cor- 
respond to the respective projections onto the Co (Fe) 
atoms (j — 1) and Pt atoms (j — 2) muffin-tin spheres. 

This subdivision yields e^cX ( 6 mca) an< ^ 6 mca- Bc- 
cause the total e^ CA contains also the contribution from 
the interstitial region between the muffin-tin spheres, we 
cannot expect that the sum e^cA + e MCA ^ s equal to 
e MCA- I n the following we call e M 'c A the "direct con- 
tribution" of the Pt atoms to the magnetic anisotropy 
energy, because it is directly assigned to the Pt sites. 
Within the framework of the rigid band model we can 



determine these quantities as a function of the band fill- 
ing Z + AZ, where Z denotes the band filling for the 
real system. The change in the band filling leads to 
the shifts in the Fermi levels appearing in cq. (1) as: 

p(l,2) , A p (2,l) 

Z + AZ = J E ^ +AEf (e) de. The results of 

this analysis are presented in Figs. 4,5. 
(b) In addition to the calculation with SOC for all atoms, 
we performed calculations where we switched off the SOC 
term for the Co(Fe) atoms or for the Pt atoms, respec- 
tively. Thereby we changed the hybridization between 
the electronic states at various sites and hence the eigen- 
values for the two orientations. Switching off the SOC 
term at say the Pt atoms does not mean that there is no 
longer the CjJfcA contribution because there is still a dif- 
ference in the Pt-projected density of states for the two 
magnetization directions due to the hybridization of the 
Pt-localized orbitals with the Co-localized orbitals. And, 
vice versa, when switching on the very large SOC term 
at the Pt atoms the hybridization of the Co-localized or- 
bitals with the Pt-localized orbitals changes drastically 
and this has a large effect on e^Q A . In the following 
we call this the "indirect contribution" of the Pt atoms 
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FIG. 5: The calculated magnetic anisotropy energy cmae = 
e(6 = 0°,(p = 0°) - e(0 = 90°, <j> = 0°) for the Co (top) 
and the Fe (bottom) wire; total contribution (+), Co (Fe) 
contribution (x), Pt contribution □, free-standing wire (■). 



to the magnetic anisotropy energy. The results of such 
analysis are presented in Figs. 6,7. 




FIG. 6: The calculated magnetic anisotropy energy £mae = 
e(6 = 0°,(p = 22.85°) - e{6 = 0°,</> = -67.15°) for the Co 
wire (top) and e M AE = e(fi = 0°, <j> = 80.40°) - e(6 = 0°, = 
—9.60°) for the Fe wire (bottom) projected to the Co (Fe) 
sites with SOC at all sites (+) or SOC just at the Co (Fe) 
sites (x), and projected to the Pt sites with SOC at all sites 
(□) or just at the Pt sites (■). 



III. DISCUSSION 

By looking at Fig. 3 it becomes immediately clear that 
there is a big influence of the Pt substrate. Whereas for 
the free-standing wire the easy axis is in the wire direc- 
tion, it is perpendicular to the wire direction for the Pt- 
supported systems. Furthermore, there is a big difference 
between the Pt-supported Co and Fe wires concerning 
the orientation of the easy axis in the plane perpendicu- 
lar to the wire, see Fig. 2. Whereas for the Fe wire the 
easy axis is nearly perpendicular to the surface normal, it 
is inclined by an angle of 22.85° from the surface normal 
towards the step edge in the case of the Co wire. Re- 
member that for the chains of Co atoms at the step edge 
of Pt(997) an angle of 43° was found experimentally^. 
The difference in the energy between the easy and the 
hard axis in the plane perpendicular to the wire is about 
1.3 meV per Co or per Fe atom, whereas the experimen- 
tal value is about 2 meV per Co atomi. In the following 



we want to distinguish between the "direct contribution" 
and the "indirect contribution" of the Pt atoms on the 
magnetic anisotropy by means of the methods (a) and 
(b) of Section 2. 

From Figs. 4,5 it becomes obvious that there is a strong 

ft Pt 

"direct contribution" e M ' CA to the total anisotropy. Over 
a wide range of the Z values, the contributions of the 
Co(Fe) atoms and the Pt atoms are large but of opposite 
signs. For the Co wires the "direct contributions" from 
Co and Pt atoms nearly cancel each other, and the sum of 
the two contributions is close to the total anisotropy en- 
ergy, indicating a small influence of the interstitial region. 
For the Fe wires the two direct contributions are much 
smaller than for the Co wire but the total anisotropy 
energy is of the same order of magnitude. The sum of 
e MCA ano - e MCA i s drastically different form e^ CA , indi- 
cating an unexpected large contribution of the interstitial 
region. 

The plots in Figs. 2,3 clearly demonstrate that the 
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FIG. 7: The calculated magnetic anisotropy energy cmae = 
e(6 = O°,0 = 0°) - e(6 = 90°, = 0°) for the Co (top) and 
the Fe wire (bottom) projected to the Co (Fe) sites with SOC 
at all sites (+) or SOC just at the Co (Fe) sites (x), and 
projected to the Pt sites with SOC at all sites (□) or just at 
the Pt sites (■). 



magnetic anisotropy of the Pt-supported Fe wires is 
strongly dominated by the SOC at the Pt atoms. In 



the absence of this coupling, only a small magnetic 
anisotropy originating from the SOC at the Fe sites re- 
mains. In contrast, to obtain the magnetic anisotropy of 
the Pt-supported Co wires, we have to take into account 
the SOC for both the Pt and the Co atoms. Switching 
on the SOC term at the Pt atoms modifies both e^CA 

, ft,Fe \ j ft.Pt rpv. A-c 4.- c ft. Co / ft.Fc \ 

( e MCA/ and e MCA- The modification of e MCA (e MCA ) 
represents the "indirect contributions" of Pt atoms (see 
(b) of section 2). Figs. 6,7 represent the contribution 
e MCA (top) an d e MCA (bottom) together with the re- 

ft Pt 

spective contributions e M ' CA on the one hand calculated 
with SOC switched on for both Co (Fe) and Pt atoms 
(corresponding to the curves shown in Figs. 4,5), and on 
the other hand calculated with SOC on the respective 
other atoms switched off. Obviously there is a big "indi- 
rect contribution" of the Pt atoms, i.e., a big change of 
e MCA ( e MCA.) w hen the SOC on Pt is switched off. The 
"indirect contribution" thereby is of the same order of 
magnitude as the "direct contribution". 

The Pt substrate has an influence also on other mag- 
netic properties of the wires. For instance, we have con- 
sidered a freestanding Fe biwire and Fe biwires at the 
step of a Pt substrate, both for a ferromagnetic and for 
an antiferromagnetic alignment of the magnetic moments 
on the two adjacent monatomic wires (within each wire 
the alignment was ferromagnetic) . The energy difference 
between these two configurations was about —20 mRy 
per Fe atom for the free-standing biwire and only about 
—7.7 mRy per Fe atom for the Pt supported biwire. For 
comparison, the difference in energy between the fer- 
romagnetic and the antiferromagnetic configuration in 
bec Fe at the equilibrium lattice constant of Pt is about 
—25 mRy/atom. The magnetic moment of the Fe atom 
for th Pt-supported biwire was very similar to the one of 
the Pt-supported monatomic wire. 

We are indebted to G. Bihlmayer and A. B. Shick for 
helpful discussions and for providing their data and re- 
sults. 
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